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Abstract 

A formal symplectic structure on R x M is constructed for the unsteady flow 
of an incompressible viscous fluid on a three dimensional domain M. The 
evolution equation for the helicity density is expressed via the divergence of 
the associated Liouville vector field that generates symplectic dilation. For an 
inviscid fluid this equation reduces to a conservation law. As an application 
the symplectic dilation is used to generate Hamiltonian automorphisms of 
the symplectic structure which are then related to the symmetries of the 
velocity field. 
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The helicity which is first discovered in has been recognized to be 
an important ingradient of the problem of relationship between invariants 
of fluid motion and the topological structure of the vorticity field 
For three-dimensional flows its ergodic and topological interpretations were 
introduced and investigated in f5|-[|l(J. It has also been studied in the context 



of Noether theorems |~T|-|T7|. Kinematical aspects of helicity invariants in 



connection with the particle relabelling symmetries were discussed in |lj 

In this work, we shall show that there is also a dynamical content of the 
helicity density in the sense that the information contained in the Eulerian 
dynamical equations can be represented in the framework of symplectic ge- 
ometry by a current vector field governing the dynamics of helicity. More 
precisely, starting from the Navier-Stokes equations of incompressible fluids 
we shall construct helicity four-vector whose divergence will define the time- 
evolution of helicity density. The dynamical properties of the fluid, such as 
viscosity, are implicit in this vector field. The evolution equation for the 
helicity density reduces to a conservation law for inviscid Euler flows. For 
fluid dynamical content of this work we shall refer to [[| and the necessary 
mathematical background can be found in ||19||-||21||. 

The Navier-Stokes equation for a viscous incompressible fluid in a bounded 
domain M C -R 3 is 

— + v- Vv = -Vp+z/V 2 v (1) 
at 

where v is the divergence-free velocity field tangent to the boundary of M, p 
is the pressure per unit density and v is the kinematic viscosity. The identity 
v ■ Vv = |V|v| 2 — v x (V x v) can be used to bring the equation into 
the form 

- v x (V x v) = z/V 2 v - V(p + -v 2 ) (2) 
and in terms of the vorticity field w = V x v this gives 

— V x (v x w) = z/V 2 w . (3) 

dt v ' w 

For a fluid with a potential <p and velocity field v the densities 

H = • V x v , H w = ■ V x w , q — w • Vy? = w((p) (4) 
will be called helicity, vortical helicity and potential vorticity, respectively. 
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Proposition 1 For a velocity field satisfying Eqs. and for q ^ 2v7i w the 
two-form 

Vl v = — (Vy? + v x w — z/V x w) • dx A dt + w • (dx A dx) (5) 

zs symplectic on I x M where I is an open interval in R. Moreover, it is 
exact, fl u = d9 with the Liouville (or canonical) one-form 

9 = -( V + p+^v 2 )dt + v -dx (6) 

which is independent of the viscosity v. 

Proof: VL v is closed by Eq. (f|) and the divergence-free property of the 
vorticity field. The non- degeneracy follows from the recognition that the 
density in the symplectic volume Q u A £l u /2 is the function q — 2v7i w which 
is assumed to be non-zero. The exactness can be verified using Eq. (Q). • 

For an arbitrary smooth function / of (t, x) the unique Hamiltonian vector 
field Xf defined by the symplectic two-form (|) via i(Xf)(Q u ) = —df is given 
by 

X f = —L-— [-^o(f)(± + v) + fw + ((Vy? — z/V x w) x V/) ■ V] . (7) 

Here, d/dt denotes the convective derivative dt + v • V which, viewed as a 
vector field on / x M, is not Hamiltonian. In fact, with the notation v = v- V, 
one can check that the one-form 

i(p t + v)(n v ) = (Vy? - z/V x w) ■ (dx - vdt) (8) 

is not closed and hence dt + v is not even locally Hamiltonian. 

Next proposition describes invariantly the connection between the sym- 
plectic structure (|5]) and the the helicity density. 

Proposition 2 The identity 

d(9 Afi„)-O„Afl,E0 (9) 

gives the equation 

™ + V • (Wv + i(p - ^v 2 )w) = U ■ V 2 w - vH w (10) 
for the evolution of helicity density. 
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Proof: We have fl u A Q u = — 2(q — 2u7i w )dx A dy A dz A dt and we compute 

V • [<pw + v x (Vf - iA7 X w)] = 2q - 2vH w - w ■ V 2 w (11) 

for the derivative of certain terms in the expression 

6„ A ft„ = 2H dx A dy A dz - 

[{<p + p- r« 2 )w + 2Hv + v x (V(f - zA7 x w)] ■ dx A dx A dt (12) 

for the three-form. Putting them together in the identity (§) we obtain Eq. 
fllUp. Upon integration, the term vw ■ V 2 w/2 in Eq. ( jTD| ) gives the integral 
of —v7i w and one obtains the usual expression for the time change of total 
helicity as given in, for example, Ref. 0. • 

Note that the helicity flux in Eq. (0) is independent of the function 
ip which we have introduced by hand to make the symplectic form non- 
degenerate. 

Using the invariant description (§) of the evolution of helicity density, 
we shall introduce a current vector J v and show that it is an infinitesimal 
symplectic dilation of Q v . J v will be defined as the one-dimensional kernel 
of the three-form 9 A Q u . Since the symplectic two-form is nondegenerate, it 
can be obtained as the unique solution of 

i(J v )(n v An v /2)=0ASl v , (13) 

that is, as the dual of the three-form 9 A Vt„ with respect to the symplectic 
volume. We find 

J v = ^rr [ZH(d t + v) + ( V + p- \v 2 )w + v x ( - vV x w) ■ V] (14) 

q — 2urL w Z 

as the expression for the helicity current. 

Proposition 3 J v is a vector field of divergence 2 with respect to the sym- 
plectic volume and it is an infinitesimal symplectic dilation for Q u . The 
evolution of helicity density 7i can be described by the identity 

div Q M -2 = 0. (15) 
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Proof: The exterior derivative of Eq. ( |T3| ) gives 



di(j v )(n v An v /2) = Cj v (Q v /\n v /2) = div nv {j v )n v /\n v /2 (16) 

= d(0Afi„) = Sl v ASl v (17) 

where we used the identity Cj = i(J) o d + d o i{ J) in the first equality and 
the second equality is the definition of the divergence. We see that J u is a 
vector field whose divergence is 2. From the last equality, we conclude that 
the equation (jl5f) is equivalent to Eq. flTOl ) describing the evolution of helicity 
density. J v is the unique vector field satisfying 

i{j v ){n v ) = 9 (is) 

and it follows from this that J v fulfills the condition 

c Jv (n v ) = di{j v ){n v ) = dd = n v (19) 



of being an infinitesimal symplectic dilation for Vt v [22]. J v is also called to 
be the Liouville vector field of Q u plf . • 

We observed that the local existence of a Hamiltonian function for d t + v 
is being prevented by the viscosity term p3[. Moreover, the viscosity term 



causes the helicity not to be conserved. We shall now show that, for the case 
of inviscid incompressible fluids described by the Euler equation, namely Eq. 
([I]) with v = 0, dt + v is Hamiltonian and that the helicity density 7i is 
conserved. To this end, we assume that the scalar field ip is advected by the 
fluid motion 

^+v-V^ = (20) 
and that the potential vorticity q ^ 0. 

Proposition 4 j^J Let v and ip satisfy Eq. ffi) with v — and Eq. 
respectively. Then, the suspended velocity field d t + v on I x M and q~ lr w are 
Hamiltonian vector fields for the exact symplectic two-form 

Q = -(Vy? + v x w) • dx A dt + w ■ (dx A dx) = d6 (21) 

with the Hamiltonian functions ip and t, respectively. The evolution equa- 
tion ftH\) reduces to the conservation law in divergence form for the helicity 
density. 
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Proof: Using Eq. ( pOD dt + v can be written in Hamiltonian form i(d t + 
v)(p,o) = —dip. More generally, the Hamiltonian vector field with the sym- 
plectic two-form (^Tj) for an arbitrary function / on I x M is given by 

*/ = ^""(/X^ + v ) + f w + x V /) • V l ( 22 ) 

which clearly reduces to d t + u for / = ip and to q~ l w for f = t. The 
conservation of helicity density is obvious. • 

For the inviscid flow of the Euler equation the helicity current takes the 
form 

Jo = ~[2H(d t + v) + {<p + p- \v 2 )w + vxVfV] (23) 
q z 

while the canonical one-form remains to be the same. That means, the 
difference between the dynamics of fluid motion with v = and v ^ is 
contained in the helicity current. Thus, the dynamical content of the helicity 
is encoded in its current and this, in turn, is connected with the symplectic 
structure on/xM which was constructed as a consequence of the Eulerian 
dynamical equations. 

The realization of dynamics of fluid motion in the symplectic framework 
is useful in the study of the geometry of the motion on M and of the hyper- 
surfaces in I x M defined by the time-dependent Lagrangian invariants, that 
is, the invariants of the velocity field. The present framework also provides 
geometric tools for the investigation of scaling properties of the fluid motion 
because the action by the Lie derivative of helicity current on tensorial ob- 



jects corresponds to infinitesimal scaling transformations [2TJ. Leaving the 
discussions of these issues elsewhere, we shall conclude this work with an ap- 
plication to the symmetry structure of the velocity field which is also related 



to the results presented in [24 



Proposition 5 Let Xf be a Hamiltonian vector field for Q u . Then, the 
vector fields {Cj u ) k {Xf)^ k = 0,1,2,... are infinitesimal Hamiltonian auto- 
morphisms ofQ u . 

Proof: The symplectic two-form is invariant under the flows of Hamiltonian 
vector fields because Cx f (^u) = di(Xf)(Q u ) = d 2 f = where we used the 
identity Cx = ° d + d o i(X) for the Lie derivative, dfl u = and the 
Hamilton's equations i(Xf)(Q v ) = —df. It then follows from the identity 

£[j v jc f ] = £j» ° £x f - £x t o C Ju (24) 
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evaluated on Q u that [J V) Xf\ also leaves fl u invariant. Replacing Xf with 



[J U} Xf] in Eq. (24) we see that one can generate an infinite hierarchy of 
invariants of the symplectic two-form Q u . To see that these are Hamiltonian 
vector fields we compute 

i([j v ,x f ])(n v ) = CjMXf)m)-i(x f )(jC Jv (n„)) (25) 

= -d(J v (f)-f) (26) 



where we used Eq. (|T9|). Thus, [J v ,Xf] is Hamiltonian with the function 
Jv{f) — f- By induction one can find similarly that (Cj v ) 2 (Xf) is Hamiltonian 
with (J u ) 2 (f) — 2J u (f) + f and so on. Interchanging J u and Xf in the identity 
(H) we also obtain i{X f ){9) = J v {f). • 

In particular, we let v = 0, / = t so that X t = q~ x w and consider the 
infinitesimal Hamiltonian automorphisms {£,j ) k {q~ 1 w) 1 k — 0, 1, 2, ... of Qq. 



The identity ( ]2"4D evaluated on the vector field d t + v gives 

£>[j Q ,q-^w){dt + v) = -C q -i w ([J , d t + v]) (27) 

where the vector field [Jo, dt + v] is, by proposition (5), Hamiltonian with the 
function Jo(^p) — (p = p — v 2 /2. By the Lie algebra isomorphism [Xf,X g ] = 
Xy^} defined by the symplectic structure fl , the right hand side of Eq. (p7|) 
is a Hamiltonian vector field with the function 

{t,p-\v 2 } = -w{p-\ 2 ) . (28) 



On the level surfaces defined by the constant values of the function ( |2"%D we 
have [[Jo, q^w], d t + v ] — In fact, if we restrict to the constant values of 
the function p — v 2 /2 the hierarchy of Hamiltonian automorphisms of Qq can 
be identified as the infinitesimal symmetries of the velocity field. This can 
be seen by replacing q~ l w with [J , q~ l w] in Eq. (^7\j. We thus proved that 

Proposition 6 For the Euler flow, the hierarchy of infinitesimal Hamilto- 
nian automorphisms (Cj ) k (q~ l w), k = 0,1,2,... of Q generate infinites- 
imal time- dependent symmetries of the velocity field on the level surfaces 
p — v 2 /2 = constant. 



As a matter of fact, the function p — v 2 /2 is related, in Ref. [p5|| , to the 



invariance under particle relabelling symmetries of the Lagrangian density of 
the variational formulation of the Euler equation. 
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